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Abstract. 

We obtain the classification of integrable equations of the form ut = u xxx + f(t,x,u,u x ,U xx ) using 
the formal symmetry method of Mikhailov et al [A.V.Mikhailov, A.B.Shabat and V.V.Sokolov, in What 
is Integrability edited by V.E. Zakharov (Springer- Verlag, Berlin 1991)]. We show that all such equations 
can be transformed to an integrable equation of the form vt — v xxx + f(v,v x ,v xx ) using transformations 

t, u, v, u x , v x ) = 0, and the u xx dependence can be eliminated except for two equations. 
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1. Introduction. 

The existence of "formal symmetries" has been proposed as an integrability test by Mikhailov et al [1]. 
A formal symmetry is a truncated expansion of the recursion operator in inverse powers of D = ^ , hence the 
equations admitting a formal symmetry are candidates for those equations admitting a recursion operator. 
The classifications of evolution equations ut = F(x,u,u x ,u xx ,u xxx ) (except for a class essentially nonlinear 
in u xxx ), and u t = h(x, t)u 3 + f(x, t, u, u x ) are given respectively in [1] and [2]. In this paper we obtain the 
classification of evolution equations u t — u xxx + f(t, x, u, u x , u xx ) using the existence of canonical conserved 
densities pi, i — 0, ...,3. We show that all such equations can be transformed to an equation without 
explicit x and t dependence using first order differential substitutions. The list of integrable equations 
given in [1] is further simplified by the availibility of time dependent transformations. Furthermore, the 
u xx dependence is eliminated, but two of the transformed equations are nonlocal. Thus, up to first order 
differential substitutions, there are only three equations in the class considered. 

2. Formal symmetries, recursion operators and differential substitutions. 

Let Ut — F[u] be an evolution equation where F is a differential polynomial. A symmetry a of this 
equation is a differential polynomial which satisfies the linearized equation <j t = F'o where F' is the Frechet 
derivative of F. The recursion operator R is a linear operator that sends symmetries to symmetries, i.e. Ra 
is a symmetry whenever a is a symmetry [3]. In general the recursion operator R is an integro-differential 
operator, and satisfies (R t + [R,F'])<t = 0. If the symmetries are dense in the solutions of this equation, 
then R can be defined as the solution of the operator equation 



In 1 + 1 dimensions the integral terms in the recursion operator can be expanded in inverse powers of D = ^ , 
and Eq.(2.1) can be solved recursively, leading to linear first order differential equations for the coefficients of 
R. The solvability of these equations in the class of local functions gives certain conserved density conditions, 
(pi)t — Drji. These conserved density conditions in turn give nonlinear partial differential equations for F, 
whose solutions lead to the classification. 

We consider the classification of evolution equations 



R t + [R,F'] =0. 



(2.1) 



u t = u 3 + f(x,t,u,u 1 ,u 2 ). 



(2.2) 
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In the following, subscripts denote differentiation with respect to x. The classification is obtained using the 
existence of conserved densities p , . . . ,p 3 , and we have chosen the coordinate transformations so as to put 
the equations in the form given in [1]. 

We shall use the following general result for the elimination of u 2 . Consider the evolution equation 
u t = u m + f(u m -i, . . . ,u,x,t). If this equation is integrable then p = df/du m -\ is a conserved density, 
hence J pt is a local function. The dynamical variable s defined by s = u m -\e~ -f p " , satisfies an equation 
of the form s t = s rn + g(s m -2 ■ ■ ■ , s, u, . . . u m -\) i.e. s m _i is eliminated. If po is itself a total derivative, 
then s is a local function. Otherwise, one has to introduce a new dynamical variable v = J po, and apply 
the transformation to the new equation. This transformation of course introduces non-locality, and the 
the crucial point is to ensure that the transformed equation is local. We shall also use the coordinate 
transformation x — > x + h(t)t to eliminate the term h'{t)u\. Other differential substitutions will be described 
as needed. 

3. Classification results. 

For equations of the form (2.2), the first four nontrivial conserved density conditions are given below. 
We shall obtain a classification based on the existence of these conserved densities. 



Po =—, (3.1a) 
1 (df\ 2 df 



3 \ u 2 J oui 

P2 =^l + -l2L2L + \ (3i c ) 

du 27 \ U2 J 3 du2 dm 3 
P3 =f7i, (3. Id) 

where (p;) t = Drji, i = 0, 1. 

The method is to compute the time derivative of the piS and integrate by parts until a term that 
involves the highest derivative non-linearly appears. Then the coefficient of this highest derivative has to be 
set to zero, and we obtain partial differential equations for /. 

As a first step, it can easily be seen that / has to be a third order polynomial in u 2 , hence, we obtain 
the u 2 dependence of /, and the evolution equation has the form 

ut = u 3 + A(x, t, u, ui)«2 + B(x, t, u, ui)u 2 + C(x, t, u, m). (3.2) 

Then, we obtain the following nonlinear differential equation for the u\ dependence of A, 

d 2 A 8A, 16,, 



For 4^0, the solution is given by 

A 4 Z ' du\ °" (3 ' 4) 

Hence 

Z = a(x, t, u)u\ + b(x, i, u)u\ + c(x, t, u). (3-5) 
In Sections 4, 5 and 6 we shall study respectively the cases where A = B = 0, A = 0, B ^ and A ^ 0. 

4. Classification for the case A = B = 0. 

We will show that the equations in this class consist of the KdV, mKdV, their potential forms and 
the Callegero-Degasperis-Fokas equation. Actually, the classification of integrable equations of the form 
u t — h(x, t)u xxx + f(t, x, u, u x ) has been obtained in [2], hence we will only give an outline of the derivations 
for completeness. The use of the coordinate transformations (4.2) given below is crucial in obtaining the 
classification. Using the conserved density conditions we first obtain 

ut — «3 — —k(t) 2 ul + 0.2(2;, t, u)u\ + a\(x, t, u)u\ + oq(x, t, u). (4-1) 

The first branching is given by k{da 2 /du) = 0. The coordinate transformations leaving the form of (4.1) 
invariant are u = a(t)u + 3(x, t), and 



---f fc - 2 — a 2 + 01 



u t =u 3 - X^Ml + Wl 
2 a z 



" r^x j -z 1 

x — fe +-a 2 
2^ a 



Ui 



2 a 2 a 



- —(u-/3)+0t- Pxxx + \^k 2 + ^a 2 - 3 x a x + a a. (4.2) 
a 2 a 2 a 

For fc ^ 0, we set fe = 1, 02 = 0, and using the conserved density conditions we obtain the Callegero- 
Degasperis-Fokas equation 

«t = «3 - \u\ + ui[cie 2 " + c 2 e- 2u }. (4.3) 

For k = 0, we first obtain u t = it 3 + a 2 (x,t,u)u 2 + a\(x,t,u)u\ + ao(x,t,u), where a 2 — b {x,t), a\ — 
c 2 (x, t)u 2 + Ci(x, t)u + c (x, t), a = ds(x, t)u 3 + d 2 (x, t)u 2 + d\(x, t)u + d (x, t). The branching is given by 
the condition bo c 2 = 0. For b ^ we have the potential KdV equation, 

u t = uz + ul, (4.4) 

For b = 0, using the conserved density conditions and appropriate coordinate transformations we obtain 
the KdV 

u t =u z + uui (4.5) 
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and mKdV equations 

u t — u 3 + u 2 m. (4.6) 

It is well known that all the equations in this class are related via Miura transformations, hence up to first 
order differential substitutions there is a single equation in the class u t — u 3 + f(x, t, u, u x ). 

5. Classification for the case A = 0, B ^ 0. 

We will show that all equations in this class are linearizable. The conserved density conditions imply 
that the evolution equation has to be in the form 

u t — u 3 + bi(x, t, u)u\u 2 + b {x, t, u)u 2 + c 3 (x, t, u)u\ + c 2 (x, t, u)u\ + c\(x, t, u)u\ + co(x, t, u). (5.1) 

Then the transformation u — > <j>(u) gives b\ — > ^ — 3^77 and we can set b\ = 0. Then the allowable 
transformations are of the form u — ► a(x, t)u + (3(x, t). The conserved density conditions give dc 3 /du = 0, 
and we set C3 constant by choosing a. Then we have the branching &0C3 = 0. If b a = 0, (5.1) is independent 
of U2, hence we set c 3 = 0. Then we have 

bo =b 02 {x,t)u 2 + b m (x,t)u + b 0Q (x,t), 

c 2 =c 21 (x,t)u + c 20 (x,t), 
1 2 

c i =3 & 02 m4 + ^boibo 2 u 3 + ci 2 (x,t)u 2 + cu(x,t)u + c w (x,t). (5.2) 
We have two subcases depending on whether b 02 = or 6 2 7^ 0. 

The case 602 7^ 0: In this case we use coordinate transformations u — > a(x, t)u + j3(x, t) to set 602 = 3 and 
&01 = 0. Then we obtain the following equation. 

u t — u 3 + 3u 2 u 2 + 9uuj + 3u 4 ui + Cio(ar, t)u\ + ^(cio)xit. (5.3) 

This equation can be linearized by the following sequence of Miura and coordinate transformations. Recall 
that the transformation (2.5) eliminates u 2 but in this case J po is not a local function. Thus we first use 
the potentiation U\ — ► u 2 to obtain 

3 Uo O / \ 7/ \ 

u t = u 3 - h iu 2 ui +u 1 J r u\cw{x, t) + d(t). 

4 u\ 

Then using the transformation u — > e 2u we obtain 

3 u 2 

u t = U3 - -— + cw{x,t)u\ + 2d(t)u. 
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Finally the transformation u — ► we obtain the linear equation 

ut — u 3 + ci (x, t)u\ + ^u((cio)a: + 2d). (5.4) 

The case b 02 = 0: In this case we use point transformations to set 601 = 3 and fe o = 0. Then the conserved 
density conditions give 

u t — u 3 + 3uu 2 + 3t*i + 3u 2 ui + uic w (x, t) + (c w ) x u + c Q0 (x, t). (5.5) 

The potentiation u\ — > u gives 

u t — u 3 + 3uiu 2 + u\ + u\Ciq{x, t) + J c o(x, t) dx. 

Then the point transformation u — > e u gives the linear equation 

ut = us + U1C10 +u c adx. (5.6) 



For completeness, we give the interacting soliton equation of the KdV equation that we would obtain if 
we had taken c 2 = instead of bi — in (5.1) 



U2«i 3 u? 3 

u t = u 3 -3 ^o^+o mu i- V 5 - 7 ) 

u 2 u z 2 



6. Classification for the case A^O. 

When A^O, the coordinate transformation u — > </>(£, i, u) results in 

A -> (2aui + U 6 - 20 K a)[aui + ui{cf) u b - 2^ x a) + (<^c- U </> X 6 + ^a)] \ (6.1) 

Thus for a = 0, these coordinate transformations can be used to take b = 1, c = 0, similarly for a ^ 0, wc 
can take a = 1 and 6 = 0. The cases c = and c ^ has to be considered separetely in the last subclass. As 
a result A has the forms studied in Sections 6.1-3 below. In all cases the allowable cordinate transformations 
should satisfy <p x = 0. We will use the transformation 



v = J eH P0 u 2 (6.2) 



to eliminate u 2 term. 



6.1 Classification for the case A = — ju^ 1 . 

We will show that all equations in this class are either linearizable, or transformable to the mKdV 
equation. We first obtain 

1 /2 

B = an(x, t, u)u 1 / + ai 2 (x, t, u)u\ + a\s(x, t, u). 
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Then C can be obtained by solving a fourth order o.d.e., and if ai 3 is nonzero, it has u\ln{ul) and U\ln{u\) 
terms. The first branching is given by the condition aua\s = 0. If we assume an = 0, ai3 7^ 0, the conserved 
density conditions give a i3 = 0. Hence we assume an 7^ 0. 

The case an 7^ (linearizable equation): We will show that the equation with an 7^ is linearizable. 
We have 

C =a i(x,t,u)u^ 2 + ao2(x, t,u)u\ + a 03 (x, t, u)ui + a oi (x,t,u) 
'2 , , 1 



+ 



3 , 2 5/2 

u i + ganai2Ui ■ 



The conserved density conditions give (004)2; — 0, and we set ao4 = by a coordinate transformation. As a 
result for an 7^ 0, we obtain the following equations. 

(an)t —{an)xxx + (au) x ao3 + 2( a 03) 

2 11 

(ai 2 )t = - g(an)xxan + g(aii) x - ^^03, 

1 1 2 2 

(aii)w =ganai2, (012)2; = — g a n> ( a 03)« = ^(aii^aii, 

a i =2(an) x , a 02 = 0, a 04 = 0. (6.3) 

We make the differential substitution s = ip(x, t, u, ui) which gives an equation where si, term is not present. 
Elimination of the coefficient of s 2 s\ requires that 



ip = ?pi(x, t, u)v ul + tpo(x, t, u). 

Then s 2 s 2 and s 2 are eliminated by choosing tpi and 1^0 as solutions of 

3(V>i)u ~ 012^1 = 0, 3(^1)^ + 011^0 = 0. (6.4) 

The resulting equation for s is independent of u by virtue of the equations (6.3) and (6.4) and it can 
be transformed to equation in the previous section by a linear coordinate transformation, and then it is 
linearized. 

The case an = (equations linearizable or transformable to mKdV): When an = 0, we have 
(012)2; = 0, and we can make ai 2 = by a point transformation, hemce B = 0. Then we are allowed to make 
transformations of the form u — > a(t)u + (i{t). The conserved density conditions give two equations in this 
class, the first one being a constant coefficient equation. 

u t = u 3 - x u\ + a iu^ 2 + a 02 ul + a 03 Mi + a 04 , (6.5a) 
ut = 113 - ju^u 2 , + a 03 (x 1 t)ui. (6.56) 
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We will show that both equations are transformable to an equation independent of u 2 - The conserved 
density p — — |^ is trivial, the transformation (6.2) reduces toti = 2-s/ui, and we have 

u t - u 3 + > v t - v 3 . 

4 ui 

Thus we obtain respectively the following equations 

1 3 

v t =v 3 + vi [-a 02 v 2 + -a iv + a 03 ] , (6.6a) 
v t =v 3 + a 03 wi + ^{a 03 ) x v. (6.66) 

The second equation is linear and the first one is transformable to the mKdV equation. Therefore, all 
integrable equations with A = — ju^ 1 are either linearizable, or transformable to the mKdV equation. 

6.2 Classification for the case A = — . 

In this case, using the conserved density conditions we obtain 

3 Wo 1 

u t — u 3 - -— + a i(t,u) h a 02 {t, u)u\ + a 03 ui + a ±(t, u), (6.7) 

where the coefficients have to satisfy the following equations. 

(aoi)t = - {a 01 ) u a 0i + 2(a 4) u aoi, (6.8a) 

(d02)t = — ( a 02)ti a 04 — {o-0i)uuu ~ 2(do4) u a 02 > (6.86) 

=(aoi)5u + 10(aoi)«tm«o2 + 15(aoi) uu (ao2)« 

+ aoi)«[9(ao2) u « + 16ag 2 ] + 2a i[(a 2) U u« + 8(002)1^02]- (6.8c) 

The equations (6.8a-c) form a consistent system, hence (6.7) is an integrable equation. The conserved density 
Po = — I ^ is trivial and the transformation v = In Ui gives 

Zu\ 1 , , . 

u t ~u 3 + - >v t -v 3 + -v 1 . (6.9) 

Z U\ 6 

We will show that, generically (6.7) cannot be transformed to a local equation independent of zt 2 . The 
coordinate transformation u — > (p(u, t) gives 



1 

1 r , 3 



O02 "^T4 [-</»«u«0u + ^uu + <t>u a 02\, 

9 U z 

ao4 — >0t + ^«ao4- (6.10) 



Using these coordinate transformations we take ao4 = 0, which implies that a i and a 2 are independent of 
t, hence the evolution equation in independent of x and t. If we set (<Zoi)u = 0, then 



(ao2)uuu + 8(ao2) u a 2 = 0, (6.11) 
hence the integrable equation in this class is 

3 u 2 1 

iH — u 3 - -— + o i \-a 2(u)ul (6-12) 

2 Mi Ml 

where ao2 satisfies (6.11). It can be seen that when ao2 is independent of u, it is transformed to the 
Callegero-Degasperis-Fokas equation, 

v t = v 3 - ^vf + Vl [-a 01 e- 2v + 3a 02 e 2v ]. 

On the other hand, if a 2 is not constant it can not be transformed to a local equation. 

6.3 Classification for the case A — — + c(x, t, u)] _1 tii. 

In this case, using the conserved density conditions we first solve B, as a function depending on the 
derivatives of c. Then the compatibility of the differential equations for C implies that c x = 0, and we can 
set c = 1 by choosing <f> u . It then follows that B = 0. As a result we obtain 

C = a 01 {t){u\ + 1) 3/2 + ao 2 {u,t)u 1 {u 2 + 1) + a 03 {t)ut + a 04 (t). 

We can set ao4 = by choosing <j> t , and ao3 = by a linear change in x. Then we obtain the following 
equation 

t*t = «s - 2 U i( M ? + ^r lu l + a oi( u i + 1) 3/2 + ao2(«)«i(«i + 1), (6.13) 

where 

{at)2)uuu + 8(a 02 )uao2 = 0, (a 02 )uaoi = 0. (6.14) 

Thus for aoi ^ we have an equation independent of u. 

The conserved density po = ~ 3(1 + iti) _1 ui?i2 is trivial and using the transformation v = J(l + 



if) 1 / 2 u 2 = In | ui + yj\ + u\ |, we obtain 



3 Ml 2 1 o 

u t ~u 3 + - 2 u 2 ~^v t -v 3 + -uf. 



Then if a 2 is constant (6.13) is transformed to 



1 3 

v t = V3 - + ^ v i[( a oi + ao2)e 2v + (a 02 - a Q1 )e^ 2v - 2a 02 ]. 



However for a i = 0, a 2 need not be constant and (6.13) is not transformable to a local equation. 
7. Conclusion. 

We have shown that all evolution equations of the form (2.2) can be transformed to an equation inde- 
pendent of x and t. The surprising fact is that the classification of equations with explicit t dependence is 
considerably simplified by the availibility of time dependent transformations. 

It is known that integrable equations of the form u t = 113 + f(x,t,u,u x ) are all related to the KdV 
equation via Miura transformations. The equations with u xx dependence that arc not transformable to one 
of these equations are 

Zul 1 , . » 

Ut =u 3 - h aoi h a 2(u)u 1 

2 u\ u\ 
3 

u t =u 3 - -jUi(u\ + 1) _1 U2 + a 02 (u)ui(uj + 1), 

where 

(0-02) uuu + 8(ao2) u &o2 = 0. 
An equivalent form of the first equation is 

ut = u 3 - - — + — [4u 6 + kiu + k 2 ]. 

2 Ml U\ 

It has been checked that this equation do not have fifth and seventh order local symmetries. This result 
suggest the following possibilities, (i) The formal symmetry condition will not hold at a later stage, but 
the formal symmetry condition has been checked at later stages and no inconsistency has been observed, 
(ii) The formal symmetry condition will hold at all orders, but it will not be possible to write the recursion 
operator in closed form. It has been shown that such a situation occurs for a system of evolution equations 
[4]. Hence these two equations need further investigation. 
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